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A  summary  is  presented  of  our  research  progress  in  three  areas:  (1)  micrometer- 
size  droplets;  ( 2 )  three-dimensional  scalar  mapping;  and  (*3)  three-dimensional  velocity 
mapping. 

Nonlinear  optical  interactions  in  a  droplet  occur  at  remarkably  low  input  intensity 
levels  because  the  droplet  acts  as  a  lens  to  concentrate  the  input  radiation  at  a  loca¬ 
tion  just  within  the  shadow  face  and  as  an  optical  cavity  to  provide  feedback  for  the 
internally  generated  radiation.  The  following  nonlinear  optical  effects  have  been 
observed  in  single  droplets:  lasing;  (2)  stimulated  Raman  scattering  up  to  the  14th 

order;  X3)  coherent  Raman  gain  due  to  the  presence  of  another  input  wave;  (4)  phase- 
modulation  broadening  of  the  elastically  scattered  and  stimulated  Raman  scattered  radia¬ 
tion;  (5)  the  delay  time  in  generating  the  multiorder  stimulated  Raman  scattering;  and 
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19.  ABSTRACT  (Continued): 

(6)  th«-\f fective  Q-factor  of  the  droplet  cavity  based  on  a  lifetime  measurement 
of  the  radiation  trapped  within  the  droplet. 

Laser-induced  breakdown  XLIB)>within  the  droplet  occurs  when  the  rising 
portion  of  the  input  laser  pulse  causes  multiphoton  ionization,  which  is  followed 
by  cascade  multiplication.  The  resultant  plasma  within  the  droplet  transforms 
a  nominally  transparent  droplet  into  an  absorping  droplet,  and  the  remaining  por¬ 
tion  of  the  input  laser  pulse  heats  the  droplet.  Plasma  is  ejected  from  the  drop¬ 
let,  first  from  the  shadow  face  and  then  from  the  illuminated  face.  Once  the 
plasma  has  been  quenched,  the  droplet  undergoes  explosive  vaporization.'.  We  have 
measured  the  following  properties  of  LIB  and  explosive  vaporization:  (1)  the 
location  of  LIB  initiation;  (2)  the  propagation  velocities  of  the  ejected  plasma; 
(3)  the  time-averaged  electron  density  along  a  line;  (4)  the  time-averaged  atomic 
temperature  along  a  line;  and  (5)  the  shape  of  the  droplet  undergoing  explosive 
vaporization. 

Three-dimensional  scalar  mapping  was  made  possible  by  sweeping  the  laser 
illumination  sheet  and  by  recording  scattered  images  on  a  high-speed  framing 
camera.,N  Such  digitally  recorded  two-dimensional  images  from  many  sheets,  all 
recorded  in  a  time  short  compared  to  fluid  motion,  enabled  us  to  construct  iso- 
cdncentration  surfaces  in  three  dimensions.  Three-dimensional  scalar  mapping 
measurements  were  made  from  aerosol-laden  cold  flows  (via  elastic  scattering) 
and  turbulent  flames  (via  molecular  specific  fluorescence  and  Rayleigh  scattering) . 

Three-dimensional  velocity  mapping  was  achieved  by  using  a  multiple-pulse 
holographic  technique.  The  turbulent  flow  was  seeded  with  many  microballoons 
that  were  coated  with  silver  to  increase  the  elastic  scattering.  An  in-line 
holography  configuration  was  determined  to  produce  the  clearest  holographic 
images  on  a  high  resolution  sensitive  film.  In  order  to  determine  the  direction 
and  velocity  of  the  seeded  particles  in  the  flow  field,  the  illuminating  beams 
were  temporally  coded.  A  digital  image  analysis  technique  was  used  to  determine 
the  direction  and  velocity  magnitude  of  the  flow  within  a  three-dimensional  sample 
volume . 
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INTRODUCTION 


During  the  three  gears  of  AFOSR  support,  significant  progress  has  been  made  in  the 
following  research  areas: 

Micrometer-Size  Droplets.  We  have  investigated  high  energy  and  high  intensity 
laser  beam  interactions  with  single  liquid  droplets  with  radius,  a,  much  larger  than  the 
incident  laser  wavelength,  X0,  i.e.,  droplets  with  large  size  parameters  X  =  27ta/X0.  Our 
efforts  involved  nonlinear  optical  effects,  such  as  lasing  and  four-wave  mixing  pro¬ 
cesses,  and  laser-induced  breakdown  (LIB)  and  explosive  vaporization  effects  that  occur 
during  and  after  a  laser  pulse. 

Three-Dimensional  Scalar  Mapping.  We  have  developed  diagnostic  techniques  that 
can  illuminate  the  sample  volume  by  a  series  of  rapidly  scanned  thin  sheets  and  can 
record  data  from  many  points  on  each  sheet  during  a  time  short  compared  to  the  flow 
time  scales. 

Three-Dimensional  Velocity  Mapping.  Our  preliminary  results  on  three-dimensional 
velocity  mapping  in  aerosol  seeded  turbulent  flows  were  achieved  by  developing  a  multi¬ 
ple-pulse  volume  hologram  technique  that  uses  silver  coated  glass  balloons  as  markers 
and  an  image  processing  formalism  to  specify  the  locations  of  each  marker  in  the  recon¬ 
structed  images. 

RESEARCH  ACCOMPLISHMENTS 

A  brief  summary  of  our  results  in  each  of  these  research  areas  follows.  Specific 
details  of  our  accomplishments  regarding  droplets  and  scalar  mapping  can  be  found  in  the 
publications  resulting  from  the  research  (see  list  on  page  16).  All  these  papers  have 
been  submitted  to  AFOSR  in  both  preprint  and  reprint  form. 
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The  intent  of  this  phase  of  the  research  was  to  develop  laser  diagnostic  techniques 
applicable  for  two-phase  flow.  The  first  objective  was  to  study  bright  wavelength- 
shifted  markers  that  can  be  used  to  trace  out  the  flow  fields  and  be  spectrally  isolated 
from  the  unadvoidable  elastic  scattering  from  the  container  walls  and  turbulent  medium. 
The  second  objective  was  to  investigate  nonlinear  spectroscopic  approaches  that  can 
identify  the  species  within  the  droplet  (as  opposed  to  the  gas  surrounding  the  droplets). 
While  we  were  pursuing  these  two  objectives,  several  unexpected  topics  related  to 
nonlinear  optical  effects  in  droplets  warranted  investigation  by  us. 


When  a  plane  wave  is  incident  on  a  droplet  with  a  large  X,  the  droplet  can  be 
envisioned  as  a  lens  that  concentrates  the  radiation  in  three  places:  (0  outside  the 
shadow  face,  (2)  just  inside  the  shadow  face,  and  (3)  just  inside  the  illuminated  face. 
While  these  three  locations  can  be  qualitatively  predicted  by  geometric  optics,  the  exact 
amount  of  intensity  enhancement  in  the  locations  requires  detailed  Lorenz-Mie  calcula¬ 
tions.  Complete  knowledge  of  these  high  intensity  regions  is  important  to  all  nonlinear 
optical  and  LIB  studies.  Our  collaboration  with  Prof.  Peter  W.  Barber's  group  at  Clarkson 
University  resulted  in  computer-intensive  calculations  of  the  internal-  and  near-field 
intensity  distribution  of  a  large  X  sphere  and  in  an  experimental  mapping  of  these  three 
locations  of  intensity  concentration  by  a  fluorescence  technique.  See  publication  *1. 

Another  unique  feature  of  a  large  X  droplet  is  that  its  spherical  liquid-air  interface 
is  capable  of  trapping  some  of  the  internally  generated  radiation  (e.g.,  fluorescence  or 
spontaneous  Raman)  when  the  wavelength  is  equal  to  one  of  the  numerous  wavelengths 
which  correspond  to  morphology-dependent  resonances  (MDR's).  With  geometries  optics, 
it  is  known  that  internal  rays  which  strike  the  liquid-air  interface  with  angles  larger 
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than  the  critical  angle  experience  quasi-total  internal  reflection.  Using  a  physical  optics 
description,  those  internal  rays  which  maintain  the  same  phase  front  after  completing 
one  trip  around  the  entire  droplet  circumference  are  trapped  within  the  spherical  inter¬ 


face.  The  droplet  for  specific  wavelengths  acts  as  an  optical  cavity  to  provide  feedback 
for  the  internally  generated  radiation.  Because  of  the  optical  feedback,  lasing  and 
stimulated  Raman  scattering  (SRS)  have  been  observed  from  liquid  droplets  with 
astonishingly  low  pump  intensity  thresholds.  Color  photographs  show  that  the  radiation 
of  lasing  droplets  and  of  droplets  undergoing  SRS  is  confined  just  inside  the  liquid-air 
interface  where  the  optical  feedback  is  high.  See  publications  *2  and  *3  for  color 
photographs. 

The  lasing  spectra  were  noted  to  consist  of  sharp  peaks  within  the  broad  fluores¬ 
cence  wavelength  region.  These  lasing  peaks  are  separated  by  nearly  equal  wavelengths, 
in  accordance  with  the  MDR's  of  a  sphere  with  large  X  and  with  the  index  of  refraction 
m  >  1.  Small  changes  in  the  lasing  peaks  have  been  used  as  indicators  of  small  radius 
changes  resulting  from  the  evaporation  of  each  droplet  within  the  linear  stream  of 
flowing  droplets.  In  addition,  utilizing  the  small  amount  of  wavelength  oscillation  in 
these  peaks,  we  have  deduced  the  amount  of  small  shape  oscillation.  From  the  frequency 
and  the  damping  rate  of  the  oscillation  amplitudes,  we  can  extract  the  dynamic  surface 
tension  and  the  bulk  viscosity  of  the  droplet.  See  publication  *4  for  a  summary  of  this 
diagnostic  technique. 

The  SRS  spectra  can  be  used  to  provide  information  on  the  molecular  species  con¬ 
tained  within  the  droplet  and,  to  a  lesser  extent,  information  on  the  species  concentra¬ 
tion.  The  energy  loss  indicated  by  the  wavelength  of  the  SRS  peaks  (usually  measured  in 
wave  number  shirts  from  the  input  laser  radiation)  is  equal  to  the  molecular  vibrational 
energy.  Thus,  the  wave  number  locations  of  the  SRS  peaks  are  fingerprints  of  the  molec¬ 
ular  species  contained  within  the  droplet.  The  intensity  of  the  SRS  peaks  has  been  sho'«n 
to  monotonically  increase  with  species  concentration.  However,  fluctuations  in  the  SRS 


peak  intensity  are  large  and,  therefore,  the  SRS  intensity  is  not  an  accurate  measure  of 
species  concentration.  For  more  details  on  the  use  of  SRS  for  species  identification, 
see  publication  *5. 

In  addition  to  lasing  and  SRS,  three  additional  nonlinear  optical  effects  were 
observed  in  large  X  and  transparent  droplets.  One  such  nonlinear  effect  is  multiorder 
Raman  scattering  which  can  be  readily  seen  from  liquid  droplets.  The  second-order  SRS 
can  be  pumped  by  the  internal  fields  of  the  incident  radiation  and  the  first-order  SRS. 

The  third-order  SRS  can  be  pumped  by  the  internal  field  of  the  second-order  SRS.  With 
such  successive  pumping  by  the  first-order  SRS  process,  up  to  the  Mth-order  SRS  has 
been  observed.  See  publication  jr6. 

The  second  nonlinear  optical  effect  observed  from  large  X  and  transparent  droplets 
is  phase-modulation  broadening  of  the  elastically  scattered  radiation  and  of  the  SRS. 
Because  the  internal  radiation  at  X0  and  at  the  multiorder  SRS  wavelengths  is  so  intense, 
the  refractive  index  of  the  droplet  is  altered,  i.e.,  m  changes  from  m0  at  low  intensity  to 
m0  +  n2l(t)  at  high  intensity,  where  n2  is  commonly  referred  to  as  the  intensity-depen¬ 
dent  index  of  refraction.  Since  the  refractive  index  is  proportional  to  the  time-varying 
intensity  l(t),  the  internal  waves  circumventing  the  droplet-air  interface  experience  a 
time-varying  phase  shift  and,  hence,  a  frequency  modulation.  Using  a  time-averaged 
detection  system,  the  effect  of  frequency  modulation  is  observed  as  a  wavelength  broad¬ 
ening  of  the  scattered  radiation  originally  centered  at  X0  and/or  at  the  various  multi¬ 
order  SRS  wavelengths.  See  publication  *7. 

Coherent  Raman  gain  is  the  third  nonlinear  optical  effect  observed  from  large  X  and 
transparent  droplets  which  are  irradiated  by  two  input  beams  with  wavelengths  at  X0  and 
X|.  The  second  beam  provides  additional  Raman  gain  and  parametric  signals  at  the  first- 
order  Stokes  of  the  first  beam.  The  combined  effect  of  the  additional  gain  and  the  para¬ 
metric  signals  due  to  the  presence  cf  the  second  beam  is  to  lower  the  SRS  threshold  cf 


the  first  beam.  See  publication  ff8. 


Our  research  capabilities  during  the  last  year  were  greatly  improved  by  the 
acquisition  of  several  items  of  capital  equipment.  In  particular,  a  streak  camera  enabled 

l 

j  us  to  determine  the  time  delay  in  building  up  the  SRS  in  single  droplets  and  the  decay 

. 

“  time  of  the  SRS  radiation  trapped  within  the  droplet.  The  first-order  SRS  build  up  time 

is  between  5  and  7  nsec  for  Q-switched  Nd:YAG  laser  pulses  (x0  =  0.532  pm,  6  nsec 
pulse  duration)  and  is  =0  50  pSec  for  mode- locked  Q-switched  Nd:YAG  laser  pulses  (X0  = 
0.532  pm,  100  psec  pulse  duration).  The  build  up  time  of  multiorder  SRS  is  much  shorter 
than  that  of  first-order  SRS.  We  have  ascribed  this  to  the  fact  that  second-order  SRS 
does  not  necessarily  start  from  spontaneous  Raman  noise  but  can  start  from  the  para¬ 
metric  signal  created  by  the  laser  and  the  first-order  SRS.  Furthermore,  the  gain  experi¬ 
enced  by  the  second-order  SRS  is  more  than  that  for  the  first-order  SRS  because  of  the 
additional  Raman  gain  provided  by  the  coherent  Raman  mixing  effect.  (See  publication 
*9.)  The  decay  time  of  the  SRS  trapped  within  the  droplet  depends  on  the  Q-factor  of 
the  mor's  of  the  droplet,  which  acts  a?  an  optical  cavity  to  provide  the  feedback  for  the 
internally  generated  Raman  radiation.  The  SRS  decay  time  is  noted  to  arise  from  two 
factors:  (1)  the  leakage  of  the  SRS  out  of  the  droplet,  commensurate  with  the  Q-factor 
of  the  MDR's  that  can  be  calculated  from  the  Lorenz-Mie  theory  for  a  perfect  sphere,  and 
(2)  the  intensity  depletion  of  the  SRS,  resulting  from  the  partial  nonlinear  optical  con¬ 
version  of  the  first-order  SRS  photons  to  second-order  SRS  photons,  which,  in  turn,  are 
partially  converted  to  third-order  SRS  photons.  By  measuring  the  decay  time  of  the  SRS 
photons,  the  effective  Q-factor  of  the  droplet  cavity  was  deduced.  See  publication  *  10. 

Several  invitations  to  international  conferences  on  nonlinear  optics  gave  us  the 
opportunity  to  review  our  findings  on  optical  effects  from  single  droplets.  See  publica¬ 
tions  *11  -  *  M. 

Laser-Induced  Breakdown  (LIB) 


In  laser  diagnostics  of  droplets  in  a  two-phase  flow,  the  upper  limit  of  the 
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laser  intensity  that  can  be  sent  into  a  samole  chamber  is  set  by  the  LIB  threshold.  The 
presence  of  fuel  droplets  greatly  reduces  the  upper  limit  of  the  laser  inensity  as  a  result 
of  the  focusing  effect  of  each  droplet  and  the  lower  breakdown  threshold  of  some  liquids 
relative  to  the  surrounding  gas.  For  this  reason,  we  decided  to  investigate  the  LIB 
characteristics  of  transparent  droplets  with  large  X. 

Ethanol  droplets  containing  efficient  fluorescent  molecules  are  known  to  reach 
the  laser  threshold  at  remarkably  low  input  pump  intensity.  We  therefore  investigated 
the  emission  properties  of  lasing  droplets  at  input  pump  intensity  levels  far  beyond  the 
lasing  thresholds.  The  following  sequence  of  emission  events  occurs  as  a  function  of 
increasing  pump  intensity:  (1)  lasing  emission  is  reached  with  low  input  intensity; 

(2)  both  SRS  and  laser  emission  are  noted  at  higher  input  intensity;  (3)  the  internal  laser 
intensity  is  sufficient  to  induce  SRS  at  still  higher  input  intensity;  and  (4)  the  LIB 
threshold  is  reached  when  the  input  pump  intensity  reaches  *1  GW/cm2,  causing  the 
appearance  of  an  intense  continuum  emission  associated  with  the  recombination  of  the 
plasma  and  the  deceleration  of  the  electrons.  Once  the  LIB  threshold  has  been  reached, 
the  plasma  emission  is  observed  from  not  only  within  the  droplet  but  also  outside  it  in 
the  region  behind  the  illuminated  face.  Whereas  the  emission  within  the  droplet  consists 
of  a  broad  continuum  that  extends  toward  the  UV  range,  the  emission  outside  the  droplet 
consists  of  a  broad  continuum  and  discrete  peak  emission,  which  is  characteristic  of 
fluorescence  from  highly  excited  atomic  and  ionized  species,  such  as  singly  ionized 
nitrogen  N(ll)  and  singly  ionized  oxygen  0(11)  in  air.  For  more  details,  see  publication 


One  of  the  key  questions  in  the  field  of  the  LIB  of  transparent  droplets  with  X  >>  ' 
is  the  location  of  the  breakdown  initiation,  i.e.,  is  the  breakdown  initiated  within  the 
droplet  shadow  face  where  the  internal  intensity  is  concentrated  or  outside  the  droplet 
shadow  face  where  the  external  intensity  is  even  higher  than  the  internal  intensity7  We 
modified  our  spectrograph  in  such  a  way  that  the  spatial  information  along  the  vertical 
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spectrograph  slit  is  preserved.  If  a  two-dimensional  optical  multichannel  vidicon  is 
placed  at  the  exit  plane  of  the  spectrograph,  this  vidicon  can  detect  both  the  spectral 


information  dispersed  along  the  horizontal  direction  and  the  spatial  information  pre¬ 
served  along  the  vertical  direction.  This  instrument  then  enables  us  to  determine  simul¬ 
taneously  the  plasma  emission  from  many  points  (both  within  and  outside  the  droplet) 
along  an  axis  parallel  to  the  laser  beam.  Our  results  indicate  that  LIB  is  initiated  out¬ 
side  the  shadow  face  for  water  droplets  surrounded  by  a  low  LIB  threshold  gas  such  as 
Ar.  However,  LIB  is  initiated  inside  the  shadow  face  for  water  droplets  (radius  <  40  pm) 
surrounded  by  air.  In  the  case  of  larger  water  droplets  (radius  >  50  pm)  surrounded  by 
air,  the  location  of  LIB  initiation  shifts  to  the  region  outside  the  shadow  face.  For 
fluorocarbon  (C8F18)  liquid,  known  for  its  high  dielectric  breakdown  strength,  the  LIB  of 
C8F |8  droplets  is  always  initiated  in  the  gas  region  outside  the  droplet  shadow  face.  For 
more  details,  see  publications  *  16  and  *  1 7. 

Another  key  question  in  field  of  the  LIB  of  transparent  droplets  relates  to  the  elec¬ 
tron  density  and  species  temperature  within  different  parts  of  the  plasma  plumes  that 
are  ejected  from  the  droplet  interface.  We  used  the  spatial  preserving  spectrograph- 
vidiccn  instrument  described  above  to  determine  the  plasma  emission  line  shape  along  a 
line  parallel  to  the  laser  beam.  The  spectral  linewidth  of  discrete  emission  peaks  from 
atomic  species  is  known  to  be  broadened  by  the  electric  fields  of  the  electrons  sur¬ 
rounding  the  atomic  species,  i.e.,  atomic  line  broadening  results  from  the  electric  fields 
of  the  electrons  through  the  first-  or  second-order  Stark  effect.  We  were,  therefore, 
able  to  extract  the  spatial  variation  of  the  electron  density  along  a  line  by  noting  the 
linewidth  broadening  of  the  atomic  emission  peak  along  a  line  that  encompasses  the 
regions  outside  the  droplet  illuminated  and  shadow  faces.  In  addition  to  a  density  varia¬ 
tion,  an  asymmetry  in  the  electron  density  is  noted  in  these  two  regions.  Based  on  the 
appearance  of  a  line  reversal  in  the  line  shape  of  the  resonance  emission  from  Na  or  Li 


region  in  front  of  the  shadow  face  is  optically  thicker  than  the  region  behind  the  illumi¬ 
nated  face.  The  intensity  ratio  of  two  discrete  wavelength  peaks  of  the  same  species 
can  be  used  as  an  indicator  of  the  species  temperature,  provided  that  local  thermody¬ 
namic  equilibrium  has  been  reached  within  the  plasma  plume.  We  measured  the  intensity 
ratio  of  the  hydrogen  Balmer  lines  and  were  able  to  extract  the  hydrogen  temperature 
along  a  line  which  encompasses  the  regions  outside  the  droplet  illuminated  and  shadow 
faces.  See  publication  *  18. 

The  spectroscopic  instrument  used  to  determine  the  location  of  tne  LIB  initiation, 
the  electron  density,  and  the  atomic  species  temperature  had  spatial  resolution  but 
lacked  temporal  resolution.  A  streak  camera  that  has  a  time  resolution  in  the  0.1  nsec 
range  was  used  to  measure  the  propagation  velocity  of  the  plasma  emission  front  by 
imaging  a  line,  which  encompasses  the  region  within  the  droplet  and  the  regions  outside 
the  droplet  illuminated  and  shadow  faces,  on  the  vertical  slit  of  the  streak  camera. 

Using  this  technique,  we  also  determined  the  propagation  velocity  of  the  plasma  emission 
for  the  following:  (1)  the  plasma  ejected  from  the  droplet  shadow  face;  (2)  the  plasma 
traveling  inside  the  droplet  from  the  shadow  face  toward  the  illuminated  face;  (3)  with 
the  laser  pulse  on,  the  plasma  ejected  from  the  illuminated  face  after  the  internal 
plasma  has  reached  the  illuminated  face;  and  (4)  with  the  laser  pulse  off,  the  external 
plasma  traveling  toward  the  laser  from  the  droplet  illuminated  face.  At  high  input 
intensites,  the  propagation  velocities  of  the  plasma  traveling  toward  the  laser  provide 
information  on  the  optical  detonation  wave.  The  response  of  our  streak  camera  is  not 
fast  enough  for  us  to  determine  the  velocity  of  the  breakdown  wave  that  is  initiated  just 
outside  the  droplet  shadow  face  during  the  rising  portion  of  the  laser  pulse.  See  publica¬ 
tion  *  1 9. 

A  transparent  droplet  is  transformed  into  an  absorbing  droplet  once  LIB  occurs.  The 
plasma  produced  by  the  LIB  process  during  the  initial  part  of  the  laser  pulse  can  absorb 
and  scatter  the  remaining  part  of  the  incident  radiation.  Once  LIB  is  initiated  within  the 


droplet  shadow  face  region,  the  absorbed  laser  energy  is  subsequently  localized  in  this 
region.  This  localized  absorption  eventually  heats  the  entire  droplet  via  the  thermal  dif¬ 
fusive  mechanism  which  takes  several  milliseconds  for  micrometer-size  droplets.  A 
framing  camera  enabled  us  to  photograph  a  droplet  every  50  nsec  with  a  frame  time  of 
10  nsec.  Using  a  back-illuminated  technique,  we  photographed  a  transparent  water 
droplet  at  various  time  delays  after  irradiation  by  a  green  high  intensity  laser  beam 
(X0  -  0.532  jim).  Such  ultrafast  photographs  reveal  that  vapor  first  emerges  from  the 
droplet  shadow  face  and  that  the  shadow  face  is  partially  consumed  by  the  vaporization 
process.  At  a  subsequent  time,  long  after  the  laser  pulse  is  shut  off,  more  and  more  of 
the  droplet  shadow  face  is  consumed  as  the  heat  travels  from  the  shadow  face  toward 
the  illuminated  face.  We  were  able  to  measure  the  ejected  vapor  velocity,  the  receding 
velocity  of  the  liquid-air  interface  starting  from  the  shadow  face  side,  and  the  velocity 
of  the  remaining  droplet  as  it  is  propelled  toward  the  laser  by  the  rapid  vaporization 
occurring  at  the  shadow  face  after  LIB  has  been  initiated.  See  publication  *20. 

A  request  to  contribute  a  paper  for  a  special  issue  of  Applied  Optics  devoted  to 
propagation  and  scattering  in  the  atmosphere  gave  us  a  chance  to  review  our  work  in  the 
field  of  LIB  in  transparent  liquid  droplets.  See  publication  *21. 

Detector  Improvement 

During  the  previous  AFOSR  contract  period,  we  needed  to  increase  the  sensitivity 
of  the  silicon  intensified  target  (SIT)  vidicon  camera.  A  summary  of  our  results  on  the 
coupling  of  an  electrostatic  image-intensif ier  on  one  end  and  a  SIT  vidicon  on  the  other 
end  of  a  coherent  fiber-optics  bundle  can  be  found  in  publication  *22.  The  image-inten- 
sifier  raised  the  lower  detectivity  limit  of  the  standard  SIT  vidicon. 


10 


MM  «.» A* A*  Itfuyjypwl 


Turbulent  flow  Is  inherently  three-dimensional  and,  although  two-dimensional  mea¬ 
surements  have  provided  valuable  insight  into  important  large-scale  flow  phenomena, 
ambiguities  remain.  Work  has  been  done  to  extend  two-dimensional  scalar  mapping  tech¬ 
niques  into  the  third  dimension.  To  solve  the  general  problem  of  three-dimensional  mea¬ 
surements  in  turbulent  flow,  data  must  be  recorded  from  many  points  within  a  volume  in 
times  on  the  order  of  10  pisec  (as  dictated  by  the  flow  time  scales).  This  requirement 
places  high  demands  on  both  the  laser  source  and  the  data  acquisition  system  used  for 
the  measurement.  During  the  research  period,  progress  has  been  made,  however,  and 
instantaneous  three-dimensional  measurements  in  turbulent  flames  have  now  been 
demonstrated. 

In  early  phases  of  the  work  on  the  extension  of  laser  imaging  techniques  from  two 
to  three  dimensions,  several  different  approaches  were  pursued.  One  of  the  easiest 
means  of  assembling  three-dimensional  information  on  flow  structures  is  to  investigate 
an  externally  forced  repeatable  flow.  By  performing  measurements  at  a  constant  time 
delay  after  the  introduction  of  a  perturbation,  a  series  of  two-dimensional  measurements 
taken  at  different  times  can  be  assembled  into  a  three-dimensional  data  set.  This 
approach  is  desirable  since  it  removes  the  requirement  of  very  high-speed  data  acquisi¬ 
tion  and  allows  signals  from  relatively  weak  scattering  mechanisms  to  be  accumulated. 
From  the  standpoint  of  fluids,  however,  there  is  no  doubt  that  the  structures  in  forced 
flows  do  not  adequately  represent  the  characteristics  of  fully  developed  turbulent  flow. 

For  this  reason,  a  continued  effort  was  made  to  develop  techniques  for  very  high¬ 
speed  imaging  of  multiple  illumination  sheets.  In  one  experiment,  a  series  of  parallel 
sheets  was  produced  by  rapidly  scanning  a  thin  illumination  sheet  produced  by  a  cavity- 
dumped  Ar*  iaser  through  an  aerosol  seeded  jet.  The  image  of  the  Lorenz-Mie  scattering 
from  the  aerosols  was  directed  onto  different  regions  of  a  doubly  intensified  vidicon 


detector  by  a  rotating  mirror.  Since  the  duration  of  pulses  from  the  cavity  dumped  Ar*  j 

laser  was  only  10  nsec,  no  blurring  of  the  images  occurred  as  a  result  of  the  rotation  of 

I 

I 

the  mirrors.  With  this  approach,  it  was  possible  to  achieve  an  effective  framing  rate  of  i 

i 

i 

30  kHz,  thus  allowing  six  different  flow  planes  to  be  recorded  in  200  psec.  < 

The  experiment  was  subsequently  improved  by  replacing  the  rotating  mirror  in  front  j 

of  the  detector  with  a  high-speed  electronic  framing  camera.  With  this  modification,  a  j 

framing  rate  of  100  kHz  was  possible,  allowing  us  to  record  16  frames  corresponding  ■ 

to  16  different  locations  of  the  sweeping  illumination  sheet.  With  our  new  detector,  the  1 

main  limitation  now  had  to  do  with  the  relatively  low  power  of  the  laser  source.  Using 
an  Ar*  laser,  the  laser  energy  in  each  exposure  time  was  <1  jij.  Thus,  Lorenz-Mie  scat-  J 

tering  was  the  only  viable  scattering  mechanism  that  could  be  used.  A  complete  discus¬ 
sion  of  this  preliminary  work  can  be  found  in  publications  *23  and  *24. 

In  our  most  recent  work,  we  have  succeeded  in  performing  a  truly  instantaneous 
three-dimensional  measurement  in  a  turbulent  flame  by  detecting  molecular  scattering. 

The  experiments  were  done  at  the  Combustion  Research  Facility  at  Sandia  National  Labo¬ 
ratories  using  a  flashlamp  pumped  dye  laser  as  the  illumination  source.  By  using  this 
laser  in  conjunction  with  the  high-speed  electronic  framing  camera,  it  was  possible  to 
record  the  fluorescence  from  12  parallel  light  sheets  intersecting  a  volume  of  the  flow 
during  a  single  1  psec  laser  pulse.  Because  of  the  high  average  power  during  the  rela¬ 
tively  long  1  jtsec  pulse,  molecular  scattering  could  be  detected.  The  framing  camera 
was  used  in  its  highest  framing  rate  mode  of  2  x  107  frames/sec.  This  work  is  described 
in  publication  *25. 

THREE-DIMENSIONAL  VELOCITY  MEASUREMENT 

Because  it  is  the  velocity  that  enters  directly  into  the  equations  of  motion  for  fluid 


flow,  its  measurement  is  of  critical  importance.  For  this  reason,  the  use  of  multiple- 
pulse  holography  to  provide  velocity  information  at  a  large  number  of  points  within  a 
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volume  was  also  investigated.  In  this  work,  a  hologram  of  an  aerosol  seeded  flow  was 
recorded  on  film  for  later  digital  analysis.  To  provide  optimum  scattering  from  the  tag¬ 
ging  particles,  aphenolic  microballoons  were  coated  with  a  thin  layer  of  silver.  This 
allowed  sufficient  scattering  to  adequately  expose  the  high  resolution,  low  sensitivity 
film  used  to  record  the  hologram.  Following  the  development  of  the  hologram,  the  real 
image  was  reconstructed  in  a  series  of  planes  digitized  using  the  high  resolution 
scanning  diode  array  camera.  Algorithms  for  efficient  processing  of  the  large  data  set 
obtained  from  the  digitized  real  images  were  also  developed.  Some  details  of  this  work 
are  given  in  the  following  sections. 


The  two  most  commonly  used  holographic  methods  differ  primarily  in  their  means 
of  reference  beam  generation.  For  in-line  holography,  a  single  beam  is  used  as  both  the 
reference  and  object  beam.  Off-axis  holography  uses  two  separate  beams  with  the  ref¬ 
erence  beam  and  the  scattered  object  beam  being  combined  at  the  recording  medium. 

Both  approaches  were  investigated.  In-line  holography  was  found  to  provide  the  best 
results  because  of  its  relative  simplicity,  moderate  film  resolution  requirements,  and  the 
efficiency  of  forward  scattering  from  the  aerosol  particles  used  to  tag  the  flow. 

Although  the  feasibility  of  double-pulse  holography  for  determining  velocity  has 
previously  been  demonstrated,  the  method  has  not  gained  widespread  acceptance.  This  is 
due  in  part  to  the  lack  of  an  automated  data  analysis  process  to  extract  velocity  data 
from  the  reconstructed  image  pairs.  One  of  the  problems  in  the  analysis  is  the  difficulty 
in  determining  which  images  correspond  to  a  given  particle.  In  addition,  there  can  be 
ambiguity  in  determining  the  direction  of  the  flow.  Both  of  these  problems  can  be  made 
more  tractable  by  using  a  large  number  of  temporally  coded  pulses  to  form  the  hologram. 
The  larger  amount  of  information  included  in  the  hologram  can  be  used  in  the  analysis 
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phase  to  resolve  the  image  matching  and  directional  ambiguities  that  limit  the  success 
of  the  double-pulse  technique. 

To  form  the  multiple-pulse  hologram,  an  argon-ion  laser  beam  was  modulated  to 
form  a  series  of  illumination  pulses  of  the  desired  pulse  duration  and  separation.  The 
modulation  was  produced  by  two  acousto-optic  modulators  that  were  controlled  by  a 
pulse  forming  circuit.  The  duration  of  each  pulse  was  selected  so  that  the  particles  in 
the  flow  would  not  move  more  than  1/10  of  their  diameter  during  the  pulse.  The  separa¬ 
tion  between  pulses  was  varied  to  remove  directional  ambiguity,  and  both  the  pulse  dura¬ 
tion  and  separation  could  be  adjusted  to  suit  the  flow  conditions.  For  our  experiments,  a 
series  of  five  pulses  was  used.  The  pulse  duration  was  8  psec  and  the  total  measure¬ 
ment  time  (i.e.,  the  time  between  the  first  and  fifth  pulse)  was  500  psec.  The  laser 
beam  was  spatially  filtered  to  provide  a  smooth  intensity  profile  and  then  passed  through 
an  air  jet  seeded  with  aerosols.  The  photographic  film  that  recorded  the  hologram  was 
located  300  mm  from  the  jet,  ensuring  that  the  far-field  requirement  of  in-line  holo¬ 
graphy  was  satisfied. 

The  aerosols  consisted  of  silver-coated  aphenolic  resin  microballoons.  Chemically 
coating  the  particles  with  a  thin  layer  of  silver  significantly  improved  the  image  con¬ 
trast  and  resulted  in  clear,  bright  images.  The  increase  in  contrast  obtained  by  chemi¬ 
cally  coating  the  microballoons  was  observed  experimentally  and  also  confirmed  by  a 
layered  sphere  calculation.  The  microballoons  had  a  density  of  0.072  g/cm3  and  a  dia¬ 
meter  range  of  5  -  124  pm.  A  90  mm-diameter  hologram  of  the  aerosol  was  recorded  on 
Kodak  technical  film  2415,  which  provided  adequate  resolution  (320  line  pairs/mm)  and 
sensitivity  (0.2  ergs/cm2)  for  the  in-line  hologram. 

Hologram  Reconstruction 


The  optical  setup  for  the  reconstruction  process  was  similar  to  the  setup  used  in 
the  formation  of  the  hologram,  except  for  the  absence  of  the  seeded  jet  and  the  addition 


of  a  600  mm  foca!  length  lens  to  collimate  the  beam.  The  collimated  beam  illuminated 
the  hologram  and  produced  a  real  image  of  the  multiply  exposed  aerosol  particles. 

Spatial  filtering  was  used  to  remove  most  of  the  undiffracted  beam  so  that  a  high  con¬ 
trast  image  of  the  aerosol  could  be  produced.  A  lens  that  focused  the  undiffracted  por¬ 
tion  of  the  beam  to  a  point  at  the  lens'  focal  plane  was  placed  behind  the  hologram  in 
order  to  accomplish  the  high-pass  spatial  filtering.  This  undiffracted  portion  was 
blocked  by  an  opaque  spot  on  a  microscope  slide.  The  diffracted  light  from  the  hologram 
formed  a  real  image  of  the  aerosol  further  from  the  spatial  filtering  lens.  The  real 
image  of  the  aerosols  fell  directly  onto  a  high  resolution  photodiode  array  camera.  Dif¬ 
ferent  planes  of  the  image  were  brought  into  focus  by  moving  the  hologram  (mounted  on  a 
stepping  motor  driven  stage)  toward  or  away  from  the  spatial  filtering  lens. 

Image  Analysis  and  Velocity  Determination 

There  are  three  steps  in  the  procedure  to  map  the  three-dimensional  velocity  field 
using  the  multiple-pulse  hologram  of  the  seeded  flow.  First,  the  location  of  each  aerosol 
image  in  the  hologram  must  be  accurately  determined.  Second,  the  seed  images  must  be 
grouped  into  a  set  that  was  generated  by  a  particular  seed  object.  Finally,  the  velocity 
magnitude  and  direction  must  be  determined  for  each  group.  There  are  several  physical 
constraints  on  the  images  that  can  be  exploited  to  make  this  complex  problem  more 
tenable. 

The  location  of  each  image  must  be  determined  in  the  first  phase  of  the  analysis. 
The  data  contained  in  the  hologram  are  recorded  by  focusing  a  set  of  two-dimensional 
planes  of  the  hologram  onto  a  high-resolution  camera  and  digitizing  the  resulting  images. 
Two  components  (x  and  y)  of  a  particle  image  location  can  be  obtained  by  computing  the 
centroid  of  an  iso-intensity  contour.  If  the  optics  used  to  reconstruct  the  hologram  are 
selected  so  that  the  depth  of  focus  of  the  real  image  is  small,  an  accurate  measurement 
of  the  z  location  can  be  made  by  determing  the  location  of  the  best  focus  of  the  point 


image.  The  best  focus  can  be  determined  by  considering  both  the  maximum  intensity  and 
minimum  image  dimension  of  the  particle.  Upon  completion  of  this  process,  a  list  of  all 
the  image  positions  and  intensities  will  be  recorded  in  an  array. 

The  next  step  in  the  image  analysis  of  the  hologram  is  to  group  the  images  into 
sets  that  can  be  assigned  to  the  same  aerosol  particle  (5  images/particle  were  used  for 
our  experiment).  As  a  first  step,  the  problem  can  be  simplified  by  using  our  knowledge 
of  the  laser  pulsing  pattern  and  the  maximum  expected  velocity  in  the  flow.  With  this 
information,  a  reduced  volume  can  be  determined  within  which  all  of  the  images 
corresponding  to  a  single  particle  must  reside.  The  list  of  all  image  locations  can  be 
searched  and  those  falling  outside  this  volume  can  be  excluded.  If  the  aerosol  seeding 
density  is  low,  there  may  be  only  5  particle  images  in  this  volume  and  the  grouping  will 
be  complete.  If  more  than  5  images  lie  within  this  volume,  an  additional  search  can  be 
performed  to  determine  which  images  most  closely  lie  along  a  straight  line  and  satisfy 
the  constraints  imposed  by  the  temporal  coding  of  the  light  pulses.  As  groups  of  images 
are  found,  they  are  eliminated  from  the  original  table  to  make  subsequent  searches  more 
efficient. 

Once  the  first  two  steps  are  complete,  it  is  relatively  straightforward  to  determine 
the  magnitude,  orientation,  and  direction  of  the  velocity.  The  orientation  can  be  deter¬ 
mined  from  a  least  squared  straight  line  fit  through  the  points  in  three  dimensions.  The 
direction  is  determined  from  the  irregular  spacing  of  the  images  (corresponding  to  the 
temporal  coding  of  the  light  pulses)  and  the  magnitude  from  the  distance  between 


extreme  points. 
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